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The DAMA/Nal set-up of the DAM A experiment has been operative during seven annual 
cycles and has investigated several rare processes. In particular, it has been realised in order 
to investigate the model independent annual modulation signature for Dark Matter particles 
in the galactic halo. With the total exposure collected in the seven annual cycles (107731 kg 
x day) a model independent evidence for the presence of a Dark Matter particle component 
in the galactic halo has been pointed out at 6.3 a C.L.. Some of the many possible corollary 
model dependent quests for the candidate particle have been presented as well. 



1 Introduction 

DAMA is an observatory for rare processes based on the development and use of various kinds 
of radiopure scintillators. It has realised several low background set-ups; the main ones are: 
i) DAMA/Nal (~ 100 kg of radiopure Nal(Tl)), which was put out of operation in July 2002 
| l | 2 | :i | 4 | 5 | () | 7 | 8 | 9 | l() | ll | 12 | l l ii) DAMA/LXe (~ 6.5 kg liquid Xenon) 14 -; iii) DAMA/R&D, devoted 
to tests on prototypes and small scale experiments'*"^ iv) the new second generation set-up 
DAM A/LIBRA (~ 250 kg; more radiopure Nal(Tl)) in operation since March 2003. 

The results obtained by DAMA/Nal, investigating over seven annual cycles (107731 kg x 
day total exposure) the presence of a Dark Matter particle component in the galactic halo by 



"also: Universita "Campus Bio- Medico" di Roma, 00155, Rome, Italy 
'also: University of Zhao Qing, Guang Dong, China 



means of the model independent WIMP annual modulation signature, have been published in 
rc fPto gether with some of the many possible corollary model dependent quests for the candidate 
particle. We invite the reader to refer to reft^ 1 for a discussion of experimental and theoretical 
arguments related to the results and to model dependent comparisons and - more in general - 
to our literature to distinguish what DAMA has presented (and its meaning) with respect with 
what sometimes is quoted by others. 

We remind that the annual modulation signature^^ is very distinctive since it requires the 
simultaneous satisfaction of all the following requirements: the rate must contain a component 
modulated according to a cosine function (1) with one year period, T, (2) and a phase, to, that 
peaks around ~ 2 nd June (3); this modulation must only be found in a well-defined low energy 
range, where WIMP induced recoils can be present (4); it must apply to those events in which 
just one detector of many actually "fires" (single-hit events), since the WIMP multi-scattering 
probability is negligible (5); the modulation amplitude in the region of maximal sensitivity is 
expected to be <7% (6). This latter rough limit would be larger in case of other possible 
scenarios such as e.g. those in refs. I -^ I IM I To mimic such a signature spurious effects or side 
reactions should be able both to account for the whole observed modulation amplitude and to 
contemporaneously satisfy all the requirements; no one has been found or suggested by anyone 
over about a decade. 

The DAMA/Nal set-up and its performances have been descri bed in refP and some other 
information on its performances and upgrading can be found in refs l^^ l Here, we just remind 
that the two PMTs, coupled through 10 cm long Tetrasil-B light guides to each Nal(Tl) crystal, 
worked in coincidence with hardware thresholds at single photoelectron level in order to assure 
high efficiency for the coincidence at few keV level^. The energy threshold of the experiment, 
2 keV, was instead determined by means of X ray sources and of keV range Compton electrons 
on the basis also of the features of the noise rejection procedure and of the efficiencies when 
lowering the number of available photoelectrons^. 



2 A brief description of the final model independent result over 7 annual cycles 

A model independent approach on the data of the seven annual cycles offers an immediate 
evidence of the presence of an annual modulation of the rate of the single-hit events in the 
lowest energy region as shown in the residuals of Fig. ^ ~~ left, where the time behaviours of 
the measured (2-4), (2-5) and (2-6) keV single-hit events residual rates are reported. The data 
favour the presence of a modulated cosine-like behaviour (A- cosuj(t — to)) at 6.3 a C.L. and 
their fit for the (2-6) keV larger statistics energy interval offers modulation amplitude equal 
to (0.0200 ± 0.0032) cpd/kg/keV, t = (140 ± 22) days and T = | = (1.00 ± 0.01) year, all 
parameters kept free in the fit. The period and phase agree with those expected in the case 
of a WIMP induced effect (T = 1 year and to roughly at ~ 152.5 t/l day of the year). The x 2 
test on the (2-6) keV residual rate disfavours the hypothesis of unmodulated behaviour giving 
a probability of 7 • 10 -4 (x 2 /d-o.f. = 71/37). The same data have also been investigated by a 
Fourier analysis as shown in Fig. right. Modulation is not observed above 6 keV^. Finally, 
a suitable statistical analysis has shown that the modulation amplitudes are statistically well 
distributed in all the crystals, in all the data taking periods and considered energy bins. More 
arguments can be found in ref. ™. A careful investigation of all the known possible sources of 
systematic and side reactions has been regularly carried out and pu blish ed at time of each data 
release while detailed quantitative discussions can be found in refs. I21H3 C _ No systematic effect 

c We take this opportunity just to remind that the experimental set-up, located deep underground, was equipped 
with a neutron shield made by Cd foils and polyethylene/paraffin moderator; moreover, a ~ 1 m of concrete almost 
completely surr ounds the installation acting as a further neutron moderator. Note that the sizeable discussion 
reported i nEEm already demonstrates that a possible modulation of the neutron flux (e.g. possibly observed by the 




Figure 1: On the left: experimental residual rate for single-hit events in the (2-4), (2-5) and (2-6) keV energy 
intervals as a function of the time over 7 annual cycles (total exposure 107731 kg x day); end of data taking July 
2002. The experimental points present the errors as vertical bars and the associated time bin width as horizontal 
bars. The superimposed curves represent the cosinusoidal function behaviours expected for a WIMP signal with a 
period equal to 1 year and phase exactly at 2 nd June; the modulation amplitudes have been obtained by best fit. 
See also ref. 2 . On the right: power spectrum of the measured (2-6) keV single-hit residuals calculated including 
also the treatment of the experimental errors and of the time binning. As it can be seen, the principal mode 
corresponds to a frequency of 2.737 • 10 -3 d _1 , that is to a period of ~ 1 year. 



or side reaction able to account for the observed modulation amplitude and to mimic a WIMP 
induced effect has been found. As a further relevant investigation, the multiple-hits events also 
collected during the DAMA/NaI-6 and 7 running periods (when each detector was equipped 
with its own Transient Digitizer with a dedicated renewed electronics) have been studied and 
analysed by using the same identical hardware and the same identical software procedures as 
for the case of the single-hit events (see Fig. The multiple-hits events class - on the contrary 
of the single-hit one - does not include events induced by WIMPs since the probability that 
a WIMP scatters off more than one detector is negligible. The fitted modulation amplitudes 
are: A = (0.0195 ± 0.0031) cpd/kg/keV and A = -(3.9 ± 7.9) • 10" 4 cpd/kg/keV for single-hit 
and multiple-hits residual rates, respectively. Thus, evidence of annual modulation is present 
in the single-hit residuals (events class to which the WIMP-induced recoils belong), while it is 
absent in the multiple-hits residual rate (event class to which only background events belong). 
Since the same identical hardware and the same identical software procedures have been used 

ICARUS coll., as reported in the ICARUS internal report TM03-01, and expected in a MonteCarlo simulation by 
H. Wulandari et al. as reported in hep-ex/0312050 ) cannot quantitatively contribute to the DAMA/Nal observed 
modulation amplitude, even if the neutron flux would be assumed to be 100 times larger than that measured 
at LNGS by various authors over more than 15 years. In addition, it cannot satisfy all the peculiarities of the 
signature mentioned above. Finally, we also remind that the contribution of solar neutrinos, whose flux is also 
expected to be modulated, is many orders of magnitude lower than the measured rate (see e.g. Astrop. Phys. 4 
(1995) 45). 
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Figure 2: Experimental residual rates over seven annual cycles for single-hit events (open circles) - class of events 
to which WIMP events belong - and over the last two annual cycles for multiple-hits events (filled triangles) - 
class of events to which WIMP events do not belong - in the (2-6) keV cumulative energy interval. They have 
been obtained by considering for each class of events the data as collected in a single annual cycle and using in 
both cases the same identical hardware and the same identical software procedures. The initial time is taken on 

August 7 th . See text. 



to analyse the two classes of events, the obtained result offers an additional strong support for 
the presence of Dark Matter particles in the galactic halo further excluding any side effect either 
from hardware or from software procedures or from background. 

In conclusion, the presence of an annual modulation in the residual rate of the single-hit 
events in the lowest energy interval (2 - 6) keV, satisfying all the features expected for a WIMP 
component in the galactic halo is supported by the data of the seven annual cycles at 6.3 a C.L.. 
This is the experimental result of DAMA/Nal. It is model independent; no other experiment 
whose result can be directly compared with this one is available so far in the field of Dark Matter 
investigation. 



3 A brief description of some of the possible corollary model dependent quests for 
a candidate 

On the basis of the obtained model independent result, corollary investigations can also be 

pursued on the nature and coupling of the WIMP candidate. This latter investigation is instead 

model dependent and - considering the large uncertainties which exist on the astrophysical, 

nuclear and particle physics assumptions and on the parameters needed in the calculations - has 

no general meaning (as it is also the case of exclusion plots, of expected recoil energy behaviours 

and of the WIMP parameters evaluated in indirect search experiments). Thus, it should be 

handled in the ' most gene ral way as we have preliminarily pointed out with time passing in 

the past^ ^^^^^^ l Some specific details can be found in ref. - and references therein, 

with devoted discussions of experimental and theoretical aspects and evaluations. Here we only 

remind that the results summarised here are not exhaustive of the many scenarios possible at 

presen t level of knowledge, including those depicted in some more recent works such as e.g. refs. 
18 19 

Fig. |3J EJ El and [7| show some of the obtained allowed regions; details and descriptions of 
the symbols are given in ref. The theoretical expectations in the purely SI coupling for the 
particular case of a neutralino candidate in the MSSM with gaugino mass unification at GUT 
scale released!^ are also shown in Fig. |S1 
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Figure 3: Case of a WIMP with mixed SI&SD interaction for the model frameworks given in ref. 2 . Coloured 
areas: example of slices (of the allowed volume) in the plane £asi vs £,°~sd for some of the possible mw and 9 
values. Inclusion of other existing uncertainties on parameters and models would further extend the regions; for 
example, the use of more favourable form factors and/or of more favourable spin factors than the ones considered 
here would move them towards lower cross sections. For details see ref. 2 . 



We remind that the allowed regions, we report herein and previously elsewhere, take into 
account the time and energy behaviours of the single-hit experimental data an d hav e been 
obtained by a maximum likelihood procedure (for a formal description see e.g. refs . IStZlHF) which 
also determines on this basis for each model the constant part of the signal and the background 
d . In particular, the likelihood function itself requires the agreement: i) of the expectations for 
the modulated part of the signal with the measured modulated behaviour for each detector and 
for each energy bin; ii) of the expectations for the unmodulated component of the signal with the 
respect to the measured differential energy distribution and - since ref. " - also with the bound 
on recoils obtained by pulse shape discrimination from the devoted DAMA/NaI-0 datM The 
latter one acts in the likelihood procedure as an experimental upper bound on the unmodulated 
component of the signal and - as a matter of fact - as an experimental lower bound on the 
estimate of the background levels. Thus, the C.L.'s, we quote for the allowed regions, already 
account for compatibility with the measured differential energy spectrum and with the measured 
upper bound on recoils. 

Specific arguments on some claimed model dependent comparisons can be found in reiPl 
They alr ead y account, as a matter of fact, also for the more recent model dependent CDMS(- 
II) claiml^ (based on a statistics of 19.4 kg • day and on a discrimination technique), where 

d We point out to the attention of the reader that the annual modulation approach gives as experimental results 
the modulation amplitudes, while the constant part of the signal, for each considered energy interval, has to be 
extracted from the measured counting rate of the single hit events, which a priori represents the sum of the signal 
and of possible residual background satisfying the trigger condition. A natural constraint on the constant part of 
the signal (and thus on the background) arises from the measured upper limits on recoils and has to be accounted 
as well, as done since ref.^1 
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Figure 4: On the left : Case of a WIMP with dominant SI interaction for the model frameworks given in ref. 2 . 
Region allowed in the plane (mw, £o~si)- The vertical dotted line represents a bound in case of a neutralino 
candidate when supersymmetric schemes based on GUT assumptions are adopted to analyse the LEP data; the 
low mass region is allowed for neutralino when other schemes are considered (see e.g. refs. 20 - 21 - 22 ) anc ] f or every 
other WIMP candidate. While the area at WIMP masses above 200 GeV is allowed only for few configurations, 
the lower one is allowed by most configurations (the colored region gathers only those above the vertical line). The 
inclusion of other existing uncertainties on parameters and models would further extend the region; for example, 
the use of more favourable SI form factor for Iodine alone would move it towards lower cross sections. On the right: 
Example of the effect induced by the inclusion of a SD component different from zero on allowed regions given in 
the plane £asi vs mw- In this example the Evans' logarithmic axisymmetric C2 halo model with vo = 170 km/s, 
Po equal to the maximum value for this model and a given set of the parameters' values (see ref. 2 ) have been 
considered. The different regions refer to different SD contributions for the particular case of 9 = 0: osd = 
pb (a), 0.02 pb (b), 0.04 pb (c), 0.05 pb (d), 0.06 pb (e), 0.08 pb (f). Analogous situation is found for the other 

model frameworks. For details see ref. 2 . 




Figure 5: Theoretical expectations of £asi versus mw in the purely SI coupling for the particular case of a 
neutralino candidate in the MSSM with gaugino mass unification at GUT scale; the curve is the same as in Fig. 

4-left. Figure taken from ref. 20 . 



DAMA /Nal is not correctly quoted and the more recent result of the 7 a nnual cycles^ is quoted 
but not accounted for. In addition, in the particular scenario of re fP, uncertainties from the 
model (from astrophysics, nuclear and particle physics) as well as some experimental ones, 
are not accounted at all and the existing interactions and scenarios to which CDMS is largely 
insensitive - on the contrary of DAMA/Nal - are ignored. 



2 10 

10 2P 




e = 2.435 



Ife.g. SI 
contribution # 
this region 
goes down 



10 



10 



200 



V 

400 

m w (GeV) 



Figure 6: Case of a WIMP with dominant SD interaction in the model frameworks given in ref. 2 . An example 
of regions allowed in the plane (mw, C°"sd) at given 9 value (9 is denned in the [0, 7r) range); here 9 = 2.435 
(Zo coupling). For the definition of the vertical line and of the coloured area see the caption of Fig. 4. Inclusion 
of other existing uncertainties on parameters and models (as discussed in ref. 2 ) would further extend the SD 
allowed regions. For example, the use of more favourable SD form factors and/or more favourable spin factors 
would move them towards lower cross sections. Values of £o\sd lower than those corresponding to this allowed 
region are possible also e.g. in case of an even small SI contribution (see ref. 2 ). For details see ref. 2 . 
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Figure 7: Case of a WIMP with preferred inelastic interaction in the given model frameworks. Examples of 
slices (coloured areas) of the allowed volumes (£<j p , S, mw) for some mw values. Inclusion of other existing 
uncertainties on parameters and models would further extend the regions; for example, the use of more favourable 
form factors and different escape velocity would move them towards lower cross sections. For details see ref. 2 . 



4 Conclusions and perspectives 



DAMA/Nal has been a pioneer experiment running at the Gran Sasso National Laboratory 
of I. N.F.N, for several years and investigating as first the WIMP annual modulation signature 
with suitable sensitivity and control of the running parameters. During seven independent 
experiments of one year each one, it has pointed out the presence of a modulation satisfying the 
many peculiarities of an effect induced by Dark Matter particles, reaching a significant evidence. 
As a corollary result, it has also pointed out the complexity of the quest for a candidate particle 
mainly because of the present poor knowledge on the many astrophysical, nuclear and particle 
physics aspects. At present after a devoted R&D effort, the second generation DAM A/LIBRA 
(a ~250 kg more radiopure Nal(Tl) set-up) has been realised and put in operation since March 
2003. Moreover, a third generation R&D toward a possible ton Nal(Tl) set-up, we proposed in 
1996 23 is in progress. 
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